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SECTION I 


INTRODUCTION 


A. DEVELOPMENT OF PILOT CONTROL 
STRATEGY IDENTIFICATION 


Pilot control strategy or piloting technique refers to the manner in 
which the human pilot perceives and processes specific information and 
manipulates the aircraft controls based on that information in order to 
perform a given task or maneuver. The principles of manual control theory 
(Ref. 1) rest on the notion that, given an aircraft and a task to perform, 
the skilled pilot is an efficient selector and processor of available and 
essential information (cues), respectively. This leads to control actions 
by which the pilot executes a task with a defined precision in a timely 
and stable manner even in the presence of disturbances, interruptions, and 
the demands of other tasks . 

It has been possible over recent years to measure how the human opera- 
tor performs a task, especially where that task is well defined. For 
example, the single-axis tracking task with one cue and one controller 
inherently restricts the basic loop structure. Further, if the controlled 
element dynamics are essentially stationary, then a spectral treatment of 
the pilot's cue-control relationship is easily obtainable based on control 
theory in the frequency domain. In this context classical describing 
functions of pilot control strategy have been widely and sucessfully ap- 
plied. Time domain parameter identification based on modern control 
theory techniques has also received substantial attention. A recent syn- 
opsis of pilot modeling approaches can be found in Ref. 2. 

The pilot identification process can be generalized in the form shown 
in Fig. 1. The objective is to structure and quantify the human operator 
using available measurements of the actual operator stimuli and re- 
sponses. This can be accomplished in many ways, however, using any number 
of combinations of the identification techniques listed in Table 1. The 


TR- 1188-2 


1 




Figure 1 . She Pilot Identification Process 
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TABLE 1 


SUMMARY OF IDENTIFICATION TECHNIQUES 


IDENTIFICATION METHODS 

Describing function analyzer 

Finite Fourier transform cross-spectral analysis 
Finite Fourier transform Input-output analysis 
Cross-correlation analysis 
Response error 
Equation error 

Sampled-data correlation (NIPIP) 

DISTURBANCE INPUTS 

Injected test inputs 

Sum of sine waves 
Frequency sweep 
Pseudo-random binary 
Random 

Existing Inputs 

Deterministic 
Random (NIPIP) 

SOLUTION CRITERIA 

Time domain 

Maximum likelihood 
Least squares (NIPIP) 

Frequency domain 

Weighted least squares 

SEARCH PROCEDURES 

Parallel-tangent 

Davidon-Fletcher-Powell-Levenberg 

Newton-Raphson 

Random 

Simplex 

Direct solution (NIPIP) 
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NIPIP approach presented in this report represents only one set of those 
techniques, but it is a powerful and relatively uncomplicated process. 

Although the computer algorithms of NIPIP are fairly simple, the real 
challenge lies in developing the assumed form of structure of the pilot 
control strategy which is to be fitted to the data. This consists of 
defining a finite difference equation which is linear with respect to 
unknown parameters. While this is partly an art at this stage, a 
practical working knowledge of manual control theory is a valuable aid in 
obtaining successful results. 


B. OBJECTIVES OP THIS STUDY 

The main objectives of the study described herein were to develop, 
implement, and demonstrate NIPIP software for NASA Dryden Flight Research 
Facility. Several examples of existing F-8 digital fly-by-wire flight 
data were treated to serve as a guide to the pilot control strategy analy- 
sis process. Additionally, applications to a future aircraft flight 
program, the AFTI/F-16, were discussed. Recommendations were then de- 
veloped regarding a future addition of interactive computer graphics to 
NIPIP. 
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SECTION II 


IDENTIFICATION OF PILOT CONTROL STRATEGY 


The purpose of this section is to provide the reader with a pragmatic 
perspective for characterizing and identifying pilot control strategy. 
This process cannot yet be done entirely automatically either through the 
use of NIPIP or any other parameter identification algorithm. Rather the 
analyst must apply his or her own skills and experience in discovering and 
quantifying the particular control strategy employed by the pilot. It is 
essential to exploit systematically all the resources at hand in this 
exercise. 


A. TASK MODELS AND PILOT MODELS 

It is necessary to distinguish the dynamics associated with the exe- 
cution of a task from those of the pilot, per se. The task model includes 
the total maneuver specification, command inputs, inanimate parts of the 
outer or task loop structure, and the overall closed-loop response of the 
pilot-vehicle combination as well as the effect of any environmental dis- 
turbances such as gusts. The pilot model is an element in that chain 
which includes the cue-control loop structure which we call piloting tech- 
nique and the perceptual features associated with the pilot. These 
distinctions are shown in Fig. 2. 

It has been found that pilot control strategy identification is us- 
ually a multistep process. It is necessary to understand first the 
overall piloting task then to hypothesize the form or structure of the 
pilot's control strategy commencing with the exterior pilot control loops 
and progressing to the interior ones. Indispensable to this progression 
is a thorough knowledge of the vehicle dynamics. This general analytical 
procedure will be demonstrated in the several examples to be presented 
shortly. 
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Figure 2. Block Diagram of Pilot-Vehicle-Task System 




B. TRACKING TASKS AND DISCRETE MANEUVERS 


There are two main kinds of piloting tasks to which NIPIP has been 
applied: tracking tasks and discrete maneuvers . In general both 

varieties can be expected to be involved in any piloting situation or 
flight phase. For example, a landing flare is itself a discrete maneuver, 
but the pilot's inner-loop regulation of pitch is a kind of tracking 
task. A change of heading is a discrete maneuver with respect to heading 
control, but there is also an inner-loop bank-angle tracking task involved 
in supporting the heading change. However the outer loop is not always a 
discrete maneuver. For example, tracking an ILS glide slope is an outer- 
loop tracking task. If flaps are selected upon intersecting the glide 
slope and initiating the descent, the flap selection could be considered 
as an "inner loop" discrete task. Clearly each case needs to be 

considered individually; nevertheless, we shall try to be as general as 
possible . 

In dealing with either tracking tasks or discrete maneuvers, it is 
useful to concentrate first on the task model. If the task model can be 
quantified adequately, then the pilot aspects, per se, can be addressed. 


C. EXAMPLES OF TRACKING TASKS 

Tracking tasks generally involve a fairly self-evident control strat- 
egy, at least on a compensatory level. There is normally a well-defined 
command input-referenced error signal and cockpit controller. For exam- 
ple, maintenance of airspeed involves a predetermined reference speed, 
V f, and a prescribed means of adjusting speed, usually either throttle 
or pitch attitude. Maintaining V rg ^ is normally a fairly long-term job, 
and it is likely that several oscillations about V re ^ would be observed 
during a sampled interval of the task execution. In contrast, a discrete 
speed change maneuver involves the step application of a new, signifi- 
cantly different value of V re ^ and the speed transient and the settling to 
a new speed . 
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Tracking tasks can be divided into explicit and implicit varieties. 
The explicit tracking task is one in which the error signal is an obvious 
presentation. Examples include gunsight pipper errors, regulation about a 
desired airspeed, heading, course, flight path angle, glide slope, or 
altitude. Some block diagrams of explicit tracking tasks with undeter- 
mined pilot describing functions are shown in Fig. 3. 

The implicit tracking tasks are ones in which an outer control loop 
supplies the command signal, an artificial construct which is itself con- 
tinually varying. Examples would include the "tracking" of pitch and roll 
in order to support such outer loop tasks as steady regulation of or dis- 
crete changes in altitude or heading. Some block diagrams of implicit 
tracking tasks with undetermined pilot describing functions are shown in 
Fig. 4. Such implicit tracking tasks will be addressed in some of the 
examples to follow. 

The pilot may use more characteristic intervals of distance (wave- 
lengths of motion) or time (natural periods of motion) to null the error 
signal in performing a tracking task than for a discrete maneuver. 
Consequently tracking tasks are sometimes characterized by lower phase 
margin (lower closed-loop damping ratio) than are discrete maneuvers. 
Both types of tasks, however, are characterized by a bandwidth 
commensurate with the task and vehicle dynamics. 


D. DEFINITION OF A DISCRETE MANEUVER 

Several examples of discrete maneuvers have been cited, such as com- 
manded changes in heading, altitude, speed, or position. In each case it 
is the "transitory" nature of a discrete maneuver which distinguishes it 
from the "steady-state" quality of a tracking task. A simple heading 
change in cruising flight is a discrete occurrence in terms of the de- 
cision to turn, initiation of the command, the turn itself, and the 
eventual settling on to the new heading. Beyond a certain point, though, 
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b. Glide Slope Displacement or Flight Director 



Note: Explicit tracking task commands are a 

clearly definable cue or instrument quantity 


Figure 3. Examples of Explicit Tracking Tasks 
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a. Height via Pitch Attitude 


Implicit 
Tracking Task 


Explicit 
Tracking Task 



b . Localizer Displacement via Bank Angle 


Explicit Implicit 

Tracking Task Tracking Task 



c. Heading or Course Angle via Bank Angle 


Explicit Implicit 

Tracking Task Tracking Task 



d. Localizer Displacement via Heading and Bank Angle 


Explicit Implicit Y 

Tracking Task Tracking Task c 



Note: Implicit tracking task commands are based on the control strategy- 

output from an exterior loop, an artificial, pilot-derived relationship 
not directly- measurable. 

Figure 4. Examples of Implicit Tracking Tasks 
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regulation of that new heading becomes equivalent to a steady-state track- 
ing task, and there may not necessarily be a precise boundary separating 
the short-term discrete task and the longer-term regulatory task. 

One useful definition for a discrete maneuver is simply the limited- 
term transition from one tracking task to another. Thus a heading change 
would be the transition from holding one heading to holding another, or a 
decelerating approach to hover would be the transition from holding a 
steady approach speed to maintaining a steady hover over the landing area. 

In the following pages the transitory aspect of the discrete manuever 
will be quantified in terms of specific mathematical models and the param- 
eters associated with those models. 


E. DISCRETE MANEUVER MODELS 

A reasonable mathematical model of the discrete maneuver can be ob- 
tained by direct transient response analysis methods. That is, a 
characteristic equation can be evaluated starting with a set of initial 
conditions and the commands appropriate to the final condition. The first 
step is to obtain a reasonable estimate of the closed-loop response 
type. One method by which this can be accomplished is through the use of 
phase plane analysis. For example, a first-order dominant mode can be 
distinguished from a second-order one depending upon the relative 
curvature of the trajectory. (Various texts can be consulted for an in- 
depth treatment of phase plane analysis, e.g.. Refs. 3 or 4.) Reference 5 
describes how pilot training manuals can be exploited to obtain estimates 
of the discrete maneuver dynamics. 


F. AN EXAMPLE OF TOOLS FOR ANALYZING 
DISCRETE MANUEVER TASK DYNAMICS 


In the course of analyzing various discrete maneuvers, it has been 
found that identification of the dominant closed-loop response mode is 
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useful. However, as the definition of a discrete maneuver implies, the 
discrete maneuver is transient . The predominant response mode may appear 
for only a fraction of its effective wavelength, period, or characteristic 
time interval. 

One technique for identifying a mode during a limited interval is to 
use a phase plane trajectory, i.e., a plot of rate versus displacement of 
a given state. The particular state to be considered is that of the dis- 
crete command. For example, in a heading change maneuver, one would 
choose to inspect heading rate plotted against heading displacement; for 
hover position, closure rate versus range; or for altitude change, 
vertical velocity versus height . Examples of such command-loop phase 
planes are presented in the following pages . 

For a second-order response, the closed-loop damping ratio, C, and 
undamped natural frequency, o) n , can be found using rigorous parameter 
identification procedures; however, even simple phase plane estimation 
methods work well. The sketch in Fig. 5 outlines one technique that has 
been found particularly useful for a variety of discrete maneuvers. Thus 
the undamped natural frequency can be extracted from the aspect ratio of 
the phase plane . A large number of landing maneuvers were so analyzed in 
Ref. 6 . 

Here we shall summarize an example of a discrete maneuver which il- 
lustrates aspects of the general closed-loop analysis technique. 

Consider the landing flare maneuver. First, based on observation of 
the closed-loop dynamics, the basic response appears to be second order — a 
damped sinusoid . Figure 6 is a sample of a landing phase-plane which 
illustrates the second-order-like behavior, at least during the latter 
portion of the trajectory. Hence a second-order transient response start- 
ing with a given height and sink rate should yield a comparable phase 
plane. If the second-order characteristic equation is assumed to be 
h + 2CfWjh + w^h = 0, then the Laplace transform can be written as 

(s + 2C f m f s + h(s) = (s + h Q + h Q (1) 
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Height, h (ft) 



Figure 6. Sample Landing Phase Plane Trajectory 
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where h Q and h Q are the Initial height and vertical velocity during the 
flare maneuver. Thus a family of general solutions could be constructed 
from the parameters Cf and and particularized using h Q and h Q . The 
appropriate command for height would presumably be zero, and this does 
appear to agree with comparisons of the above model to actual flight 
data. An example of a DC-10 landing with the matched second-order model 
parameters is shown in Fig. 7. 

For the DC-10 landing flare, it was found in Ref. 6 that a fairly 
large sample of pilots preferred a closed-loop damping ratio of about 0.7 
±0.1 and a closed-loop natural frequency of about 0.4 ± 0.1 rad/sec. It 
should be noted that a closed-loop response with these properties tends to 
provide consistently good decay of sink rate from a wide range of initial 
conditions, from off-nominal aircraft flight conditions, or from a varia- 
tion in flare maneuver aggressiveness. 

If the closed-loop response can be evaluated as shown above, then it 
may be possible to deduce something about the pilot control strategy and 
the perceptual pathways. In Ref. 6 it was shown that the combined pilot- 
vehicle system during landing has the general properties of a lag-lead 
network. Further, using ensemble landing data and knowledge of tbe air- 
craft flight path dynamics, one can deduce the use of lead-compensated 
height variation and the existence of a significant lag or decay in ad- 
dition to the airframe flight path lag. 

A general effective lag-lead pilot-vehicle form for the landing man- 
euver is shown In Fig. 8. Assumption of such a form can be based on 
knowledge of the vehicle flight path dynamics and the deduction that the 
rate feedback or its equivalent must be involved to explain the relatively 
large amount of closed-loop damping. By expanding the closed-loop char- 
acteristic equation for this feedback system, the open-loop parameters T^ 
and T-j- can be related to the closed-loop parameters 4^ and in the fol- 
lowing manner: 
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Figure 7 • Typical Landing Maneuver ] 
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Figure 8. Block Diagram of Equivalent 
Pilot-Vehicle System for Flare 
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Hence, if true lead compensation were involved in a fixed amount even for 
varying pilot gain, there should be a trend in ensemble landing data sug- 
gested by the latter equation, namely, ensemble landing data, when plotted 
in terms of versus should have a slope equal to and an inter- 
cept equal to 1/Tj as shown in Fig. 9. Such was shown to be the case in 
Ref. 6. In fact a more detailed analysis based on this concept was con- 
ducted resulting in the suggestion of lead higher than first order 
(perhaps involving vestibular as well as visual feedback) and the indica- 
tion of a substantial lag or delay beyond that of just the airframe and 
closed-loop pitch response. For actual landings, this lag was not detri- 
mental, but for simulator landings it was excessive and could be used to 
explain the tendency for hard landings. Thus this analysis procedure 
permitted an assessment of simulator fidelity and even training effective- 
ness of the simulator through direct comparison of simulator landings with 
those made in the actual aircraft. 
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Figure 9 . Regression Line Analysis Scheme 
for Ensemble Landing Data 
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SECTION III 


APPLICATION OF ANALYSIS TOOLS 


It Is assumed that the reader has available the NIPIP user's guide 
(Ref. 7) and is familiar with Section II, Background and Theory of Opera- 
tion. Therefore the description of procedures which follows will not 
require lengthy or rigorous treatment. The main goal is to provide some 
ideas for applying analysis tools either in an exploratory fashion to 
increase insight and understanding or in a more deliberate data analysis 
mode in which an accepted model form is refined and finally quantified. 

The various features of pilot control strategy that need to be con- 
sidered and which should be eventually quantified are listed in Table 2. 
In fact, the specific objective of this section is to discuss how each can 
be addressed via NIPIP. 


A. HYPOTHESIS OF GENERAL LOOP STRUCTURE 

Effective quantitative identification of pilot control strategy re- 
quires some degree of understanding of the basic flight task or 

maneuver. A written description may be available from training manuals or 
flight test reports, or an oral pilot commentary can be useful. 

References 2 and 5 illustrate how training materials can be literally 

interpreted, not only to obtain structural descriptions of the essential 
cues and feedbacks for the task but in some cases to solve directly for 
pilot control strategy gains. 

The chief feature of the general loop structure is the command loop 

state variable, the cue which forms the outermost feedback control loop. 
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TABLE 2 

QUANTIFIABLE FEATURES OF PILOTING TECHNIQUE 

Loop structure in terms of essential cues and feedbacks 

Logical switching points or criteria 

Loop gains 

Loop compensation 

Time or spatial dependence 

Sampling or discrete control strategy 

Successive organization of perception (SOP) stage 

Closed-loop pilot-vehicle response 
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In many tasks the command loop state variable Is obvious from the task 
description : 


Task 


Command Loop State Variable 


Heading change, regulation 
Altitude change, regulation 
Airspeed change, regulation 
Straight and level flight 


Heading, if 
Altitude, h 
Airspeed, V 

Heading (lateral-directional axes) 
and altitude (longitudinal axis) 


In some tasks, especially those involving outside visual reference, 
the command loop state variable is less clear. For example, in a visual 
approach the pilot may be using a complex, ill-defined geometric construct 
based on his perspective view of the airport area. Nevertheless this 
could be approximated by a simple glide-slope-like parameter for the 
purpose of quantifying control strategy. That is to say, we may not know 
the exact way in which a pilot derives visual or motion state information, 
but we can assume that the cue is essentially equivalent to the 
corresponding true state. The perceptual distortion of the true state can 
always be added to the control law if the distortion is sufficiently 
known . 

Clues to the nature of the loop structure can be obtained from the 
closed-loop task execution response. For example, active flight path 
regulation correlated with altitude suggests the presence of an altitude 
feedback. Further strong damping with respect to altitude would suggest 
either vertical velocity feedback or its equivalent . Thus candidate loop 
structure configurations can be developed by deductive reasoning based on 
manual control theory fundamentals coupled with a mathematical model of 
the aircraft and the task. This deductive approach has been described in 
Ref . 8 . 
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It is fair to point out that at some stage of loop structure hypothe- 
sizing, the analyst is likely to be faced with a level of ambiguity among 
candidates. This ambiguity may be resolvable with further data analysis. 


B. LOGICAL SWITCHING POINTS OR CRITERIA 

FOR STRUCTURE ADJUSTMENT 

It is normal to encounter changes in the basic task loop structure 
which are a function of control or state nonlinearities, the pilot switch- 
ing to other tasks, or a change in the operating environment. Such 
changes in loop structure cannot be ignored when using an identification 
scheme such as NIPIP because of the hazard in applying an invalid model 
form to identify a portion of flight data. This is, in fact, one of the 
major hurdles to creating a truly automatic pilot control strategy identi- 
fication scheme. 

Some examples of logical switching points are given in Ref. 5 for 
turns and altitude changes. If a pilot chooses to change heading more 
than, say, 30 deg, it would be common to observe a steady turn rate limit 
(or a bank angle limit) until reaching a heading sufficiently close to 
that desired. Then the pilot would roll out of the turn with a loop 

closed on heading, per se . Figure 10 summarizes the phase plane of such 
action. 

For the above maneuver the pilot's decision to turn or roll out is 
represented by a logical switch which transitions the loop structure from 
a constant bank angle command to a heading feedback as shown in Fig. 11. 
The decision or switching logic is represented by a function of heading 
error, designated as f(<f> e ). According to one widely accepted rule of 

thumb (Ref. 9), that switch in technique would occur when the heading 
error reaches one third (or one half) the steady bank angle. For example 
for a 30-deg banked turn, the pilot might begin to roll out 10 deg before 
the desired new heading. This would then be a guide to identifying the 
end of the turn maneuver. One would then apply NIPIP first to the steady 
turn with the bank angle control loop structure shown previously in 


TR- 1188-2 


24 



Turn 

Completed 


New Heading 
Commanded 


Turn 

Rate, 

♦ 


Heading Error, 
♦ ' *c 



Figure 10. Large Amplitude Heading Change Maneuver 
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Heading, i|r 


Loop Structure for 
a Constant Bank 
30 Degrees 



Fig. 3(c) and, second, to the rollout with the loop structure shown in 
Fig. 12. 

Alternatively, if the pilot chooses to change heading by regulating 
turn rate, as in the case of making a standard rate turn by observing the 
turn rate needle, the constant 30 deg bank angle command of Fig. 11 would 
not necessarily apply, especially if the speed were varying substan- 
tially. Instead one would apply NIPIP first to the steady turn, with the 
turn rate control loop structure shown in Fig. 13; then, second, to the 
rollout with the loop structure shown in Fig. 12. 


C. INTERPRETATION OF LOOP GAINS 

AND COMPENSATION 

The identification of pilot gains and compensation can be done fairly 
explicitly with NIPIP, regardless of whether the hypothetical loop struc- 
ture for the piloting technique involves nearly periodic sampling 
operations which can be reflected explicitly in the identification process 
(see Topic E, following) or nearly continuous operations which can be 
approximated by a very short sampling interval in the identification pro- 
cess. An example of speed regulation technique via throttle control 
[Fig. 3(a) herein] was identified using a sampling strategy in Ref. 10, 
whereas an example of flight director regulation via column control 
[Fig. 3(b) herein] in the same reference was identified using a continuous 
control strategy. Usually the unknowns to be solved (the c matrix in the 
user's guide, Ref. 7) can represent continuous feedback gains, or they can 
be Interpreted in terms of effective lead or lag compensation. 

In some cases it is desirable to interpret the finite difference equa- 
tion, as solved by NIPIP, in the continuous domain, because the user may 
be more familiar with forms of compensation in the continuous domain. For 
example : 


6(n) = k^6(n-l) + k 20 (n-l) 


(4) 
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Pilot 


Aircraft 



Figure 12. Loop Structure for Rolling Out of a Turn 
and Holding Heading 


Pilot 


Aircraft 



Figure 13. Steady Turning Loop Structure for 
Turn Rate Regulation 
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can be legitimately interpreted as a first-order lag if has a value 
between zero and one. If translated using an inverse z-transform, 


6(z) (l - e aT z *) = 9(z) k(l - e a ^J z ^ (5) 


Corresponds to 


6 


-a 6 — ka0 


( 6 ) 


where 



(7) 


and 



( 8 ) 


These results are not exact, however, since there is not a one-to-one 
correspondence between the continuous and discrete domains. Several other 
transformation methods could be applied with about equal accuracy. 
Table 3 shows a number of such transformations for a first-order lag. 
Similar sets could be derived for higher order continuous systems. 
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TABLE 3 


EXAMPLES OF FINITE DIFFERENCE EQUATIONS 
APPROXIMATING A FIRST-ORDER CONTINUOUS LAG 


DIFFERENTIAL EQUATION 


x = -ax + au 
LAPLACE TRANSFORM 


x(s) _ a _ 1 

u(s) s + a 1 + s/a 

FINITE DIFFERENCE EQUATION APPROXIMATIONS 


1. Direct z-Transform 


x(z) = (l - e aT )z 1 
u(z) " (1 - e -aT z -1 ) 


or x 

n 



+ (l 



u 


n-1 


2. Tustin Transform 


x(z) 

u(z) 


a(T/2) 

1 + a( T/2) 


(1 


1 - 


1 - a( T/2) 
1 + a( T/2> 



-1 

z 


or x 

n 


1 ~ a( T/2) a(T/2) 

1 + a( T/s ) n-1 1 + a(T/2) 


(u + u ) 
^ n n-1' 


3. Half-Period Advance 

x(z) , (1 - r a(T/2) ) 4. [ e - a(T/2) - e~ aT )z -1 

u(z) . -aT -1 

1 - e z 


or x 

n 


-aT . ( . -a(T/2) 

e x , + ll - e 
n-1 ' 


) u + 
' n 


-a(T/2) 



u 

n- 
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TABLE 3 (Continued) 


4. Difference Equation and Adams Second-Order Integration 


x = -ax , + au 
n n-1 n 


x = x , + (T/2)(3x - x .1 

n n-1 v n n-1' 


or x - fl - a( 3T/2) )x , - a(T/2)x 0 + a(3T/2)u - a(T/2)u . 

n v ' n— i n-Z n n— 1 


or 


x(z) 

u(z) 


a(T/2)(3 - z l ) 


1 - (l - a(3T/2))z -1 - a( T/2) z 


-2 


5. Difference Equation and Euler Integration 


x = -ax . + au 
n n-1 n 


x = x , + T x , 
n n-1 n-1 


or x = x , - aTx „ + aTu . 
n n-1 n-2 n-1 


or 


x(z) 

u(z) 


aTz 


-1 


, -1 , ^ -2 
1 - z + aTz 


6. Simultaneous Difference Equation and AdamS Second-Order Integration 


x = -ax + au 
n n n 


x 

n 


X n-! + - VJ 


or x 

n 


1 + a(T/2) a(3T/2) a(T/2) 

1 + a( 3T/2) n-1 1 + a(3T/2) u n 1 + a(3T/2) U n-1 
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TABLE 3 (Continued) 


7. Simultaneous Difference Equation and Euler Integration 


x = -ax + au 
n n n 


x = x , + Tx , 
n n-1 n-1 


or x = (1 - aT)x , + aTu , 

n n-1 n-1 


8. Simultaneous Difference Equation and Trapezoidal Integration 


x = -ax + au 
n n n 


c = x . + T/2 (x + x . ) 
n n-1 v n n-1' 


or 


* . „ + »(T /| > („ + „ ) 

" (l + a(T/2) ) "- 1 1 + a(T/2 > " "- 1 


NOTE: Same as No. 2, the Tustin Transform 


9. Fowler's Method with Half-Period Advance 


x(z) 
u( z) 


(1 - e -W2>) + ( e -(aT/2) _ e -.T u -l 

, -aT -1 
1 - e z 


or x 


-aT 
e x 


, + (1 - e- (aT/2 >)u + (e- (aI/2) - e aT )u , 

n-1 ' n ^ 'n-1 
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TABLE 3 (Concluded) 


10. Fowler's Method with z-Transform Substitution 


-aT 

x = e x 
n 


. + (l - e“ aT )u . 
n-1 v ' n-1 


11. Difference Equation and Rectangular Integration 


x = -ax , + au 
n n-1 n 


x = x , + Tx 
n n-1 n 


or x = (1 - aT)x , + aTu 

n n-1 n 


12. Modified z-Transform with Full-Period Advance 


x(z) _ (l - e aT ) 


u(z) 


, -aT -1 
1 - e z 


or x = e 
n 


-aT 


x .+(l-e a ^)u 
n- 1 v ' n 


13. Difference Equation with Lag Halved and Half-Period Advance 


x = -2ax , + 2au 
n-1 


x = e 
n 


-2aT 


x .+(l-e aT )u +( 
n-1 v ' n v 


-aT -2aT 
e - e 


)u , 

' n-1 
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D. TIME OR SPATIAL DEPENDENCE 


If a varying control strategy is likely, then at least three means of 
detection are available. These are: 


1. Sequential finite-duration time-averaging windows 

2. Sliding time-averaging window 

3. Time or spatial variation imbedded in the general 
control strategy form. 


The first of these was demonstrated in Ref. 11 for a flight director 
tracking task. The three-minute landing approach phase was simply divided 
into six 30-sec segments. Each of these was sufficiently long to permit 
reasonable convergence of a NIPIP solution. This showed a steadily in- 
creasing pilot gain well beyond the breakout height point (see Fig. 14). 

A sliding time window-averaging provides an autoregressive moving 
average solution. This offers the potential of detecting a time varying 
strategy (if it exists) in a shorter overall run length after the initial 
settling time, because it avoids the concatenation of several subsequent 
settling intervals. 

Perhaps the most satisfactory way to handle a time or spatial varia- 
tion is to imbed it directly in the NIPIP finite difference equation. 
This is feasible, though, only if a reliable model form can be assumed. 


E. SAMPLING OR DISCRETE CONTROL STRATEGY 

Often the pilot is assumed to be a continuous controller, and this is 
sometimes justifiable. It is possible, though, to discern where the pilot 
is operating more like a sampling or discrete controller and to estimate 
the parameters of that strategy, e.g., throttle regulation of airpseed in 
Ref. 10. 
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Clues to a sampling strategy appear as repetitive discrete maneuvers 
in the time domain or in the phase plane. This characteristic is likely 
to be more prominent in some outer control loops. One example given in 
the next section involves flight path control. The vertical velocity 
command is shown to be a series of steps with approximately 15 sec inter- 
vals. The strategy thus implied is that an outer loop on visual glide 
slope error is present. Further, the visual glide slope cue is being 
sampled with an approximately 15 sec period resulting in a succession of 
discrete vertical velocity commands. 


F. SUCCESSIVE ORGANIZATION OF PERCEPTION (SOP) STAGE 

It is possible to assume pilot control strategy forms for any of the 
three stages of SOP: compensatory, pursuit, or precognitive . These forms 
are summarized in Fig. 15. 

The basic compensatory level is most frequently applied, but the 
higher stages should be considered where substantial skill and control 
coordination are involved. 


G. CLOSED-LOOP PILOT- VEHICLE RESPONSE 


If the predominant closed-loop pilot-vehicle response is of second 
order (a common occurrence) , the closed-loop damping ratio, C, and un- 
damped natural frequency, u) n , can be found using the least-squares 
parameter identification procedure incorporated in NIPIP. The usefulness 
of these closed-loop parameters for interpreting pilot control strategy 
has been illustrated previously in the example of the landing flare man- 
euver at the end of Section I. Other applications are presented in 
Refs. 12 and 13. 
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a ) Initial Phase : COMPENSATORY (Single Loop) 



b) Second Phase: PURSUIT (Mull Hoop) 



c) Final Phase : PRECOGNITIVE (Open-Loop ) 

Figure 15. The Three Main Phases 
in the Successive Organization of Perception (SOP) 
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This concludes our discussion of the quantifiable features of piloting 
technique which can be addressed by application of the NIPIP. In the next 
section we shall present the results from analyses of several sets of 
flight test data acquired at the Dryden Flight Research Facility. 
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SECTION IV 


ANALYSIS CASES 


Several sets of flight data were analyzed from the handling qualities 
investigation reported in Ref. 14. Each specific analysis is presented as 
a case study in this section. The various cases demonstrated the applica- 
tion of NIPIP software and presentation of some of the pitfalls associated 
with handling flight data. 

The following cases depict several facets of pilot-vehicle analysis 
including: 


1. Closed-loop task execution dynamics 

2. Inner- and outer-loop pilot control strategy 

3. Control of flight path in vertical and lateral planes 

4. Identification of vehicle characteristics. 


Owing to problems with the quality of flight data, a number of the 
results have a practical value in terms of pointing out difficulties or 
limitations with this or any other identification process. Where flight 
data appear to be adequate, a few of the results are meaningful and reveal 
time varying pilot gains as well as the nature of the pilot compensation. 
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A. NORMAL APPROACH AND LANDING 


1. Case 1: Closed-Loop Longitudinal Task Dynamics — Approach 

As a prelude to examining pilot control strategy, the closed-loop task 
execution response is first considered. The maneuver is characterized 
both in terms of selected time histories (Fig. 16) and a phase plane plot 
of the assumed command loop state variable, vertical velocity (Fig. 17). 
No visual approach slope indication nor other visual landing aid was em- 
ployed in these approaches. 

From the time histories it appears that the approach is composed of 
several segments each having a different sink rate command and is followed 
by the flare segment. This is more readily visible in the phase plane in 
terms of the various second-order response loops. Here the flare segment 
is treated as a response to a variable sink rate command; in Ref. 6 and 
subsequently in Case 4, the flare segment is regarded as an unforced 
response from an initial velocity, sink rate, and height to a set of 
desired conditions at touchdown. Touchdown Is inferred from the abrupt 
increase in the amplitude of estimated vertical acceleration. The quanti- 
zation of height (30 ft) was too coarse to define touchdown, and sink rate 
was not recorded. 

What is not discernible in the data is the factor responsible for the 
sink rate commands, i.e., the outer-loop flight path cue. Since the pilot 
is flying a visual approach, there is no direct record of a flight path 
error signal. Thus it is possible to log only the discrete steps. 

NIPIP can be used to identify the sink rate commands as well as the 
response parameters connected with following those commands. Thus if the 
closed-loop response were considered to be second order (e.g., based on 
Fig. 17), that is, 


h 


-2C w h — a) 
a a a 


(h - hj , 


(9) 
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then NIPiP should be able to solve for u> a> and h c provided we have 

•• • •• • 

suitable data for h , h, and h. (A description of filtering used to 
obtain h , h, h is found in Appendix A.) 

The above response function was analyzed using NIPIP for each of the 
approach segments Indicated (in Figs. 16 and 17). The results are shown 
in Fig. 18. Three solutions are plotted for most of the segments. The 
solid line describes the NIPIP solution using all data which were sampled 
at 50 Hz; the broken line describes the NIPIP solution using every fifth 
data sample (sampling rate 10 Hz) from the total length of the segment; 
hence, the legend "decimated 1/5" is given for the broken line in 
Fig. 18. First the closed-loop natural frequency tends to lie in the 
range of 0.5 to 0.6 rad/sec — in basic agreement with other data previously 
analyzed (Ref. 6). This frequency can also be observed in the phase plane 
plot using the analysis technique presented earlier (i.e., the aspect 
ratio of each phase plane loop). 

The damping ratio, C a , obtained by NIPIP is very large; however, and 

we would expect values slightly less than one from Ref. 6. The damping 

ratio solutions between one and three, shown by the solid and dashed lines 

in Fig. 18, were obtained with the more exact estimate of h. The lower 

values were obtained using a less exact, more filtered estimation for 
" • •• 
h . There is the suggestion that state variable estimation and/or h basic 

flight data quality is a factor. Since the main concern of this report is 

to outline the NIPIP analysis process and not to describe techniques in 

state variable estimation, only two simple methods of estimating h were 

chosen. Without a doubt it would be possible, using more sophisticated 

filtering and estimation techniques, to reconstruct desired states in 

order to obtain more reasonable values of the second-order response 

parameters . 

Finally the NIPIP vertical velocity command solutions appear to be 
reasonable but not as accurate as one would wish in Segments 1 through 4; 
recall that the flare, Segment 5, is treated here as a response to a 
variable sink rate command for which we have no independent "apparent" 
estimate. At the same time, improved data quality could be expected to 
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improve these estimation results as well as those for closed-loop undamped 
natural frequency and damping ratio. 


2. Case 2: Longitudinal Pilot Control 

St rategy — Approach 


This case in which the task segmentation was formulated is a natural 
progression from the first. Here the role of NIPIP is to identify pilot 
control strategy characteristics for the inner (pitch attitude) and outer 
(vertical velocity) loops. A block diagram of the hypothetical control 
technique was shown previously in Fig. 4(a), provided we interpret the 
outer loop for the explicit tracking task in terms of vertical velocity, 

• • a 

h, instead of height, h, and the command as h c = h rg £. 

Figure 19 shows the primary control and command loop state variables 

for the inner and outer loops. Also shown are the four approach segments 

and one flare segment. NIPIP was used to obtain estimates of Y p ^ 

[Fig. 4(a)] and Y . for each of these finite time segments. 

Ph 

The inner-loop solutions are shown in Fig. 20, but the results are 
confounded by a fundamental data quality problem. Note that for the first 
segment (the initial pushover onto the approach path) the identified phase 
angles for Y quickly converge on about -180 deg, thus indicating an 

r Q 

established negative feedback control loop. Very quickly, however, the 
phase shifts to nearly -90 deg and remains for the duration of the ap- 
proach. Further, the amplitude of Y shows a -6 dB/octave slope (the 

P 0 

frequencies shown are separated by approximately one octave) . Hence the 
Yp^ solution appears to be an integrator — not the form expected for the 
pilot in this case. A pure gain would be considered more likely. 


The explanation for the above behavior is believed to be the coarse- 
ness in pitch attitude quantization in the flight data records. Consider 
one of the finite difference equation forms assumed for NIPIP: 


6 


e 

n 


a. 6 +b,0 . + b 9 + 

1 e . 1 n-1 2 n-2 

n-1 


c.h . + bias 
n-1 


( 10 ) 
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Note that for any interval in which the attitude and vertical velocity are 
varying within the quantization (round off) band, the last four terms on 
the right side of the equation all act as constant, say e 0 . Thus if is 
continuing to vary, 


b. 0 , + b_0 0 + ch , + bias = e 0 = constant (11) 

1 n-1 2 n-2 n-1 0 


and 


5 = a, 6 

e 1 e 

n n-1 


+ e. 


( 12 ) 


Hence for moderately slow changes in compared to the sampling rate, 


a 


n 

l 


(13) 


and 


6 g (z) (l - z » e 0 (z) 


(14) 


where e 0 is a residual non-zero quantity related to the quantization of 0 


or 


Y (z) 
p 0 


yo 

e 0 (z) 



(15) 


Thus Y would appear as an integrator with -6 dB/octave rolloff in ampli- 

r 0 

tude and -90 deg phase — the behavior observed during much of the 
approach. Therefore the NIPIP solution, except for the short periods 
during the pushover and later in the pre-flare, appears to yield a false 
pilot control strategy model, because the pitch attitude was recorded with 
too coarse a quantization in the flight test. 
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Note that in Fig. 21 where the assumed form of Yp^ is given more de- 
grees of freedom no improvement is seen. 

Similar kinds of problems carry over into the outer loop solution 

except for the fact that the false solution of Yp^ takes on a different 

character. In this case, when the attitude begins to dwell within the 

quantization bands for long periods, it appears that Y . simply reflects 

"h 

the ratio of the steady-state attitude to sink rate, i.e.. 



(16) 


e.g., for 0 = -2 deg = -0.0349 rad 

h = -75 ft/sec (17) 

and the steady-state solutions are 


and 



= 20 log yjr = -31 dB 


(18) 


< Y . 


P 


h 


-180 deg 


(19) 


both of which are displayed quite consistently in Segments 2 and 3 of 
Fig. 22 and somewhat erratically in Segments 1 and 4. 


An increased order form for Y . (Fig. 23) does little to help except 

”h 

in Segment 4 where 0 is changing significantly, and even there only the 
phase angle appears to converge to a valid solution. 


Appendix B presents the results of a more detailed investigation of 
the effects of quantization in pitch attitude on the identification of 
pitch attitude and sink rate control strategy. 
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Pilot model forms: s e ^ = a^5e n _^ + 9c n 


8~ = b 0 + b„ 0 .+ch , + bias 

c n on 1 n-1 n-1 


Figure 22 . Case 2: Outer Loop Pilot Longitudinal Control Strategy Solution 

(Lower Order Model Form), Normal Approach and Landing 
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3. Case 3: Lateral Pilot Control Strategy — Approach 


It was hoped that a heading control strategy via bank angle control, 
like that shown previously In Fig. 4(c), might permit the Identification 
of the pilot's lateral technique during at least part of the approach. 
Figure 24 shows the available time histories for the lateral control and 
command loop state variables. As indicated in Fig. 24, a gentle turning 
maneuver was performed early in the approach, ostensibly to acquire the 
final approach course, just following the pushover shown previously in 
Fig. 19. The heading (40 record, however, shows a gradual change through- 
out the entire approach. This is typical of an approach in a cross-wind 
shear and/or an approach with a varying airspeed during the final (sup- 
posedly straight) portion of which approach the pilot is tracking the 
runway centerline perspective (and extension thereof) with an outer lat- 
eral displacement regulation loop. Recorded variations in altitude and 
true airspeed suggest that this interpretation may be valid. Neither lat- 
eral displacement nor ground speed time histories were recorded, however, 
from which to identify an outer lateral displacement control strategy, 
Yp^ . Therefore no attempt could be made to identify the pilot's three- 
loop lateral technique as depicted in Fig. 4(d) through either the 
acquisition or the final (supposedly straight) portion of the approach in 
Fig. 24. 

Notwithstanding the foregoing limitation, a heading control strategy 
like that shown previously in Fig. 4(c) was assumed for the lateral axis, 
and the NIPIP was applied to the time histories for the lateral control 
and command loop state variables in Fig. 24. 


*Given an initial condition, y , one might estimate lateral 
displacement, y, from ground speed, Vg, and heading, 4, by means of the 
equation 

y = y 0 + /q v g (C) sln ♦(€) 
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Figure 2k. Case J: Time Histories for Lateral Control Strategy, 

Normal Approach and Landing 





Three segments were chosen: the early turn maneuver (0 to 25 sec), 

the final approach (25 to 70 sec), and the flare (70 to 90 sec) which 
Includes landing at 85 sec. The Inner-loop control strategy for Y Is 
shown in Fig. 25 at three frequencies, 0.1, 0.51, and 5.77 rad/sec. Note 
that In the first and second segments, the Identified phase angles for Y 
quickly converge toward -180 deg, thus indicating an established negative 
feedback control loop. Convergence of neither phase angle nor amplitude 
is sustained, however, until in the second segment the phase shifts to 
approximately -190 deg (at 0.1 rad/sec), -230 deg (at 0.51 rad/sec), and 
-270 deg (at 5.77 rad/sec); and the amplitude shows a -20 dB/decade slope 
in the frequency decade between 0.51 and 5.77 rad/sec. Similar conver- 
gence occurs in the third segment. Hence the Y solution appears to 

P <j) 

represent a first-order lag at 0.5 rad/sec at times of 33 sec and 74 sec; 
but then, at 50 sec, the phase changes to that of an integration — not the 
form expected for the pilot in this case either. First-order lag compen- 
sation for Y at 0.5 rad/sec seems low in frequency, but may be plausible 
P(j) 

since the pilot was dividing his attention between flare and line-up 
control. 

The abrupt change at 50 sec to identify an integration was not antici- 
pated here, because roll attitude quantization is 0.02 deg (one-fifth that 
of pitch attitude) , and roll attitude exhibits some variation in 
Fig. 24. At 50 sec, however, Fig. 24 shows an interval of relative quies- 
cence in the pilot's control activity and in roll attitude. In hindsight, 
it might have been preferable to change the past value of the lateral 

control stick displacement from <5 to 5„ in the difference equation 

a n-l a n-2 


6 


a 

n 


a 5 + b j + b 6 + c + bias 

1 a , 1 n-1 2 n-2 1 n-1 

n-1 


( 20 ) 


in order to verify the tendency for a^ to converge on a value which is 

approaching unity as in Case 2. Thus we offer the additional 

recommendation to replace <5 by 6 (or 5 ; ra > 2) where 

a n-l a n-2 a n-m 

quantization appears sufficiently small, yet a^ converges to unity in a 
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Figure 25. Case 3 : Inner Loop Pilot Lateral Control Strategy Solution 

Normal Approach and Landing 



solution for strategy which would not otherwise be expected to identify 
the form of an integration. One must take this recommendation cautiously, 
because there may be some other cases, such as that of a pure gain con- 
trolled element, where an integration would represent the expected 
strategy for control in the frequency range of observation. 

The outer-loop control strategy for Y is shown in Fig. 26 at three 
frequencies: 0.05, 0.2, and 2 rad/sec. Note that there is no semblance 

of convergence on a solution for Y p in the turn segment. This fact is 
not unexpected, because the loop structure of Fig. 4(c) is inappropriate 
for a turn unless a ramp function representing the heading command is 
introduced during the identification of the pilot's strategy by NIPIP. 
The heading command was taken to be zero in applying NIPIP here, however. 

During the segment labeled "final" in Fig. 26, the solution for Y 
appears to converge on a low frequency gain of 1 deg/sec with limited lead 
compensation above a frequency of 0.2 rad/sec. For reference a gain of 
3 deg/deg is plotted in order to indicate the general agreement with the 
value inferred from pilot training guides in Ref. 5. 

During the segment labeled "flare" in Fig. 26, the solution for Y p ^ 
does not appear to converge in the interval from 70 to 85 sec, at which 
time Fig. 24 indicates that landing occurred. Since the time histories in 
Fig. 24 indicate that a small lateral correction was made between 70 and 
80 sec, it is likely that a three-loop lateral control strategy including 
an outer displacement loop should be used in future attempts to identify 
the pilot's lateral approach technique with Case 3. 


4. Case 4: Flare Maneuver 

The final segment (five) of the approach was assumed to be the flare 
maneuver. The time histories of all of the variables are contained in 
earlier figures (Figs. 16 and 19). Figure 27 shows a phase plane of h and 
h and indicates another serious flight data quality problem. Notably the 
altitude quantization band is 30 ft wide. While this did not prevent 
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Apparent convergence on a low frequency gain of 1 deg/ deg 
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Figure 26 . Case 5: Outer Loop Pilot Lateral Control Strategy Solution 

Normal Approach and Landing 
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Case A: Phase Plane for Landing Flare Analysis 

Normal Approach and Landing 
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estimation of an adequate sink rate during the approach, it does present 
problems in the flare. 

The results of closed-loop task response identification by NIPIP did 
not converge to valid solutions for damping ratio, £, and natural fre- 
quency, u) n , because of the coarse quantization in recorded altitude (see 
Fig. 28). 


Inner-loop pilot control strategy results are shown in Fig. 29 and 
outer-loop results in Fig. 30. Reasonable settling in the solution occurs 
during the initial half of the maneuver. The solutions are invalid after 
77 sec, because the altitude remains invariant within the last quantiza- 
tion band whereas the estimated sink rate continues to decrease until 
touchdown at 85 sec. Both Yp^ and Yp are practically pure gains. Inner- 

loop gain, Y , is about 0.1 (-20 dB) and is comparable to those found 
P 0 

during the approach. Outer loop gain, Y , is about 0.01 deg/ft 

* h 

(-40 dB) — lower than those values found in the DC-10 study, Ref. 6 — but 
evident in the flight records. 


B. SPOT LANDING WITH A LATERAL OFFSET 

1. Case 5: Closed-Loop Longitudinal Task 

Dynamics — Approach 

In the previous cases the pilot was flying a "normal" visual approach 
and landing. In the following four cases, the pilot was to fly a spot 
landing with a lateral offset. This consisted of an approach that was 

lined up with the edge of the runway, followed by an offset maneuver 
(initiated 100 ft above ground level) to line up with the runway center- 
line, and a touchdown at the 5000 ft marker. 

As before, the closed-loop task execution response is first con- 
sidered. The maneuver is again characterized in terms of the time history 
of the assumed command loop state variable, vertical velocity (Fig. 31). 
As seen in Fig. 31, it appears that the approach is again composed of 
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Figure 28. Case Identified Parameters for Closed-Loop Flare 
Task Analysis, Normal Approach and Landing 
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Figure 29* Case Inner-Loop Pilot Longitudinal Control 
Strategy Solution, Normal Approach and Landing 
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Figure 30. Case Outer-Loop Pilot Longitudinal Control 
Strategy Solution, Normal Approach and Landing 


TR- 11 88-2 


63 



t (sec) 


Figure 31. Case 5: Time Histories for Longitudinal Closed-Loop Analysis 

Lateral Offset Approach, Spot Landing 



several segments each having a given sink rate command. These segments 
are indicated In Fig. 31. Here it is necessary to Infer touchdown from 
the abrupt Increase In the frequency of estimated vertical acceleration, 
because the quantization of height (30 ft) was too coarse to define 
touchdown, and sink rate was not recorded. 

NIPIP is again used to identify the sink rate command as well as the 
response parameters, £ a and ui a , by considering the closed-loop response to 
be second order. The second-order response function of Eq. 9 was analyzed 
using NIPIP for each of the approach segments indicated in Fig. 31. The 
results are shown in Fig. 32. In the first segment the identified solu- 
tion for the vertical velocity command, h c , does not predict the 
appropriate value; in fact, the predicted h c value was greater than zero 
beyond 7 sec. (The positive h c scale is not shown.) 

In the second segment the identified undamped natural frequency, m a , 
approaches zero during the initial 10 sec of the segment, and NIPIP esti- 
mates a subsidence and divergence thereafter. This leads to the 
estimation of an infinite value of damping ratio, C a , and zero for the 
vertical velocity command. Possible reasons for these anomalous results 
appear in the time histories of Fig. 31 between 18 and 22 sec. The pilot 
appears to have trimmed the aircraft in the descent so that pitch attitude 
is virtually constant (-2.5 < 0 < -2 deg), vertical acceleration is flue- 

•• r\ 

tuating about null (-4 < h < 4 ft/sec^), and vertical velocity is almost 
• •• • •• 
constant (-75 < h < 65 ft/sec). The estimate of h from h in this portion 

of the second segment is evidently also fluctuating about null so as to 

cause the identification of w a in Eq. 9 to approach zero. (For example, 

• • 9 •• 

u>_ must vanish in Eq. 9, if h = h = 0.) Anomalous results such as this 

a 

are typical of trimmed flight conditions where neither the pilot nor the 
turbulence is disturbing the recorded variables sufficiently to permit 
reliable estimation of the variations in the states. This further indi- 
cates that there are basic underlying limitations in using the flight 
data. Again, since the main concern of this report is to outline the 
NIPIP analysis process and not to describe techniques in state variable 

•• t 

estimation, only a single simple method of estimating h was chosen. 
(Vide Appendix B.) Provided that the trimmed flight condition is 


TR- 1188-2 


65 




disturbed sufficiently, it would be possible, using more sophisticated 
filtering and estimation techniques, to reconstruct desired states in 
order to obtain more reasonable values of the second-order response 
parameters. 

Finally in Segments 3 and 4 of Fig. 32, the overall trends in the 
identified solution are similar to those observed in Case 1. The 

identified closed-loop undamped natural frequency, o> a , is in the range of 
0.4 to 0.6 and the damping ratio, £ a , is greater than one. The vertical 
velocity command solutions appear to be reasonable for these segments. 


2. Case 6: Longitudinal Pilot Control 

Strategy — Approach 


Figure 33 shows the primary control and command loop state variables 
for the inner and outer loops. Also shown are the four approach seg- 
ments. The flare segment will be discussed in Case 8. NIPIP was used to 
obtain estimates of Y_ and Y . for each of these finite time segments. 

Pe Ph 

The inner loop solutions are shown in Fig. 34. As before, the results are 
confounded by the fundamental flight data quality problem experienced in 
Case 2, and the identified solution in each case appears to be an 
integration. 


This result Is explained using the reasoning presented in Case 2. It 
is believed to be caused by the coarse quantization in pitch attitude in 
the flight data records. Thus, as before, the NIPIP solution yields a 
false pilot control strategy model. 



For this case, additional degrees of freedom in the assumed form of 
were tried and, as before, did not improve the identified solution. 


The identification of the outer-loop pilot control strategy in Fig. 35 

also suffered from similar problems experienced earlier in Case 2. It 

appears that Y . simply reflects the ratio of the steady-state attitude to 
H h 

sink rate. Increased order for Yp^ did little to help in identifying the 
pilot's control strategy. 
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trategy. 







3. Case 7: Lateral Pilot Control 

Strategy — Approach 


As previously mentioned In Case 5, the lateral task during this ap- 
proach differed from the first by the offset maneuver. In terms of the 
NIPIP solution for lateral control strategy, only two segments were used — 
course acquisition and offset maneuver. Figure 36 reveals the recorded 
variables as functions of time and identifies the two segments in which 
NIPIP was applied. A 2 to 3 deg/sec turn was in progress during the first 
half of the course acquisition segment. The second half of the course 
acquisition segment is labeled "final" where tracking of the final ap- 
proach course commences and continues. One would hope that the heading 
control strategy via bank angle in Fig. 4(c) would yield valid solutions 
during at least the "final" half of the course acquisition. 


Lateral inner-loop control strategy solutions for Y in Fig. 37 dur- 
ing the initial turn start to converge on a first-order lag at 0.5 rad/sec 
and then shift to represent an integrator at about 15 sec. (Recall that 
-180 deg is the negative feedback reference for judging phase angle solu- 
tions for Y .) This unexpected result was also observed in the former 

P<j> 

Lateral Case 3 in Fig. 25, although formerly during the straight final 

approach and not during the initial turn. Then in Fig. 37 when the final 

course has been acquired at 20 sec, the solution for the phase angle of 

Y_ changes abruptly to represent the expected lead-lag compensation. 

P «J> 

Again at 30 sec there is a hint of another change in the amplitude of Y . 

p <t> 

The solutions for Y in the "final" half of the course acquisition 

p <j> 

(i.e., the on-course portion) were recalculated using a segment between 20 
and 45 sec. The results are shown in Fig. 38. Convergence appears to 

occur at 32 sec. The amplitude of Y is about -30 dB at low frequencies 

p 4* 

and decays at -10 dB/decade in the frequency decade from 0.51 to 
5.77 rad/sec. The phase angle of Y , which is -187 deg at 0.1 rad/sec, 

P 4> 

changes to -200 deg at 0.51 rad/sec and to -215 deg at 5.77 rad/sec. 
(Recall that the negative feedback reference angle for Y is -180 deg.) 

P ijj 

Clearly this is lag compensation, albeit slight, and not lead-lag compen- 
sation as Fig. 37 appeared to suggest. This unexpected low frequency lag 
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Figure 38. Case 7 : Inner-Loop Pilot Lateral Control Strategy 

for On-Course Segment 


(circa 0.5 rad/sec) in the estimated roll control strategy may have 
resulted from the division of attention between flare control and line-up 
control coupled with the apparent absence of disturbing variability in 
crosswind affecting lateral displacement, heading, and roll attitude. 
Consequently less frequent roll attitude corrections were needed during 
this particular on-course segment. The solutions for Y during the 

ly 

"offset maneuver" in Fig. 37 appear to converge only after 60 sec to 
represent a pure gain of -30 dB and a time delay of about one-quarter 
second, both of which are at least plausible. 


Heading control strategy solutions for Y in Fig. 39 during the 
initial turn appear to be incorrect, because the large turn rate command 
was not incorporated as a dependent variable. Even after course acquisi- 


tion occurs at 20 sec, the solutions for Y 


although changed, still 


appear to be incorrect. We would expect a pure gain on the order of 10 dB 


(3 deg/deg) for Y between 20 and 45 sec. 


The solutions for Y during 
the "offset maneuver" in Fig. 38 do not appear to converge at all. 


The solutions for Y_ in the "final" half of the course acquisition 

P* 

(i.e., the on-course portion) were recalculated using a segment between 20 
and 45 sec. The results are shown in Fig. 40, which is the counterpart of 
Fig. 38 for Y . Convergence appears to occur at 32 sec. The amplitude 

p <j> 

of Y is about +10 dB (3 deg/deg) at low frequencies (0.05 and 

0.2 rad/sec) and rises only about 2 to 4 dB in the frequency decade from 

0.2 to 2 rad/sec; the phase angle is about zero at low frequencies and 

rises only about 10 to 20 deg on the same frequency decade (0.2 to 

2 rad/sec) . Thus Y is nearly the pure gain expected with a very slight 

p <l» 

tendency for lead compensation, and it exhibits a slightly increasing gain 
with time to over +12 dB (4 deg/deg). These results appear to be plau- 
sible and compatible with the value of 3 deg/deg inferred from pilot 
training guides in Ref. 5. 
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4. Case 8: Flare Maneuver 

As shown In Figs 41 through 44, the results are similar to those ob- 
tained in the former (Case 4) Figs. 24 through 27. The large (30 ft) 
height quantization band evident in the phase plane of Fig. 41 precludes 
reliable estimates of the undamped natural frequency and damping ratio of 
the closed-loop flare in Fig. 42. The solutions failed to converge in the 
cross hatched interval in Fig. 42. The undamped natural frequency ulti- 
mately converged on a value of approximately 0.25 rad/sec, which is 
comparable to other flight results in Ref. 6; but, as in Case 4, the damp- 
ing ratio converged in the neighborhood of 0.2, an unexpectedly low value 
for such a protracted floating flare. 

The inner- and outer-loop pilot control strategy solutions are shown 
in Figs. 43 and 44, respectively, and are practically pure gains. Atti- 
tude gain, Y , is about 0.3 (-10 dB), which is triple that value found 
” 0 

previously for the flare in Fig. 26 (Case 4), perhaps because the pilot 

was to fly a spot landing in the present case. Yet the outer loop gain, 

Y , is again identified in Fig. 43 as 0.01 deg/ft (-40 dB) (cf. Fig. 27, 
"h 

Case 4) before the solution becomes invalid, because the altitude remains 
invariant within the last quantization band, whereas the vertical velocity 
remains negative until touchdown. 


C. VEHICLE IDENTIFICATION— CASES 9 AND 10 

This topic describes and demonstrates two possible techniques for 
identifying the transfer functions between the perturbed body-axis normal 
acceleration, a z , and the body-axis pitch acceleration, q, due to the 
elevator, <5. The first technique employs specific forms for transfer 
functions between the vehicle accelerations and controls with undetermined 
parameters. NIPIP can then be used to determine the unknown parameters. 
This technique is similar to the pilot identification method demonstrated 
above and in Ref. 7. Using the "short-period approximation," Ref. 15, the 
required transfer functions are given below: 
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Figure ^2. Case 8: Identified Parameters for Closed-Loop Flare 
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( 21 ) 



A 


a 

z 
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»(s + 1/T fl ) 


q s 2 + 25 W S + (to ) 2 
sp sp sp 


( 22 ) 


The z-transform equivalents of Eqs. 21 and 22 are used as inputs to NIPIP; 

A A 

the outputs of NIPIP are frequency responses of a z /6 and q / 6 . 

The problem with this first technique, however, is the tacit 
assumption that all changes in the vehicle accelerations are due to the 
specified control deflections. If the vehicle is being disturbed by gusts 
or by another control the results will be erroneous. This was proven to 
be true when this technique was applied to flight test data for the F-8 
aircraft described in Ref. 7. The transfer function technique could not 
identify a z /<5 or q/6 when the simulated pilot-vehicle system was evidently 
disturbed by random turbulence. Because this technique was unsuccessful, 
no presentation of the results is made here. The second technique, 
described below, was successful in identifying coefficients of equations 
describing the F-8 aircraft from the same flight test data; therefore we 
will present results only for the second technique. 

The second technique is to specify forms for the equations of motion 
of the vehicle and then to use NIPIP to determine the unknown 
parameters. The following relations between total body-axis normal 
acceleration, A z , body-axis pitch acceleration, Q, and the aircraft states 
are used as inputs to NIPIP, 


A 

z 



Z a + Z a 
UW WW 


Z 6 
U5 


+ Z 


UU 


(23) 
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(24) 


^2 ' V + "u, § + % 7 + V + M ira 

In Eqs. 23 and 24, U is the airspeed, a is the angle of attack, 6 is the 
elevator, and the coefficients are called "total" derivatives (Ref. 16). 
The total derivatives are related to the small perturbation stability and 
control partial derivatives* as follows: 


z 

u 

— 

2U [Z._. + Z .6 + Z TTTT W 

o L UU Uo o J UW o 

(25) 

Z 

w 

= 

Z._.U + 2Z. W 

UW o WW o 

(26) 

Z 6 

= 

vl 

(27) 

M 

u 

= 

2D o["du + Vj + Vo + Vo 

(28) 

M 

w 

= 

Vo 

(29) 

M 

q 

= 

Vo 

(30) 

M. 

w 

= 

K 

(31) 

M 6 

= 

v» 

(32) 

U 

0 

= 

V_ cos a 
I o 

(33) 



o 


W 

o 

= 

V_ sin a 
T 0 

(34) 


o 


*The partial derivatives are defined in Ref. 15. 
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The outputs of NIPIP (i.e., the "total" derivatives in Eqs. 23 and 24) 
are used as inputs to Eqs. 25 through 34. Frequency responses of a z /6 and 
q/6 can then be obtained by using approximate factors (Ref. 15) to 
calculate the gains, damping ratios, and undamped natural frequencies of 
Eqs. 21 and 22. 

The time histories of aircraft states shown in Fig. 45 were used to 
demonstrate the vehicle identification technique described above. The 
data are from a formation flight task using the NASA DFRF F-8 aircraft. 
The pilot's task was to acquire and then maintain a constant distance 
between his own and another aircraft. 

The climb and acquisition of altitude during the first 50 sec of the 
time history was not part of this task and was therefore ignored. NIPIP 
was then set up to provide estimates using a fixed window every 10 sec 
between 60 and 80 sec. Table 4 summarizes the small perturbation 
stability and control partial derivatives at 60, 70, and 80 sec. As 

stated previously, these partial derivatives were then used to compute 
transfer functions via approximate factors, which in turn were used to 

A A 

compute the frequency responses for a z /6 and q/6 shown in Figs. 46 and 47, 
respectively. The figures show a well-damped short period (Sgp = 0.5) 
with a natural frequency of about 1.5 rad/sec. The pitch numerator zero, 
1/Tg^, is about 0.5 rad/sec. 

Figures 46 and 47 also show estimates of DFRF F-8 aircraft frequency 
responses based on wind tunnel data with the SAS on. Note that the two 
sets of frequency responses are not in very good agreement. A detailed 
investigation into the reason for these differences was considered beyond 
the scope of this report. However we did compare the two different 
estimates of transient responses to the actual F-8 aircraft responses in 
the time domain. The comparison is summarized in Fig. 48. Note that the 
NIPIP estimate of q and a n due to 6 is much closer to the actual F-8 
response than the wind tunnel estimate. 
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TABLE 4 


ESTIMATED STABILITY AND CONTROL DERIVATIVES FOR THE DFRF F-8 AIRCRAFT 

(Flight 49, Form 1) 


00 


t 

(sec) 

N 

*2, 

R? 

q 

Z u 

(1/sec) 

Z w 

(1/sec) 

z 6 

(ft/sec Z /rad) 

( 1 /sec-f t) 

(1/sec-ft) 

M q 

(1/sec) 

m 6 

(l/sec z -rad) 

% 

(1/ft) 

60 

100 

0.996 

0.930 

-.0021 

-.422 

102.2 

-.0045 

-.0057 

-1.35 

-7.1 

0.0003 

70 

200 

0.996 

0.909 

-.0433 

-.376 

103.4 

-.0044 

-.00593 

-1.29 

-6.84 

0.0002 

80 

300 

0.996 

0.868 

-.0580 

-.356 

103.4 

-.0042 

-.00597 

-1.10 

-6.40 

.000015 




LEGEND: Solid line is q/5 via NIPIP using DFRF F-8 

flight data shown in Fig. 45 
Dashed line is q/5 via wind tunnel data 
with SAS on 


Figure 46. Estimated a_/s Frequency Response for the DFRF F-8 Aircraft 




TR-1188-2 



flight data shown in Fig. 45 
Dashed line is q/s via wind tunnel data 
with SAS on 


Figure 47. Estimated q/s Frequency Response for the DFRF F-8 Aircraft 
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Elevator from t = 8.5 to 13.5 sec 
for Flight k-9, Form 1 (see Fig. 45) 



Actual q response of DFRF F -8 



q response due to elevator based 
on NIPIP estimates 



q response due to elevator based 
on wind tunnel data (pitch SAS on) 



Actual a n response of DFRF F -8 



a^ response due to elevator based on 
NIPIP estimates 



a^ response due to elevator based on 
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Figure 48. Comparison of Estimated and Actual Pitch Rate 
and Specific Force Responses Due to Elevator 
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SECTION V 


RECOMMENDATIONS FOR APPLYING NIP IP 


This project was centered around the development of a user-oriented 
software package and the exercising of that software using actual flight 
data. At this point we shall present recommendations for applying NIPIP 
to other programs along with suggestions for enhancing the present soft- 
ware package . 


A. GENERAL RECOMMENDATIONS 

Regardless of the specific application, the pilot-vehicle tasks, com- 
mands, and external disturbances must be sufficient to excite the relevant 
states of the vehicle and to require pilot control activity. In applying 
NIPIP the analyst should remember that a trimmed aircraft usually reveals 
little or nothing about the pilot's control strategy, because the pilot 
is, by definition, not actively involved in the control process after the 
aircraft has been trimmed. 

The analyst of pilot control strategy should always start with 
suitable mathematical models of the task(s) and the controlled element 
before attempting to interpret the results of the NIPIP. This preparation 
not only increases the likelihood that the relevant candidates for the 
control loop structure will be exposed but also prepares the analyst with 
rational estimates for ranges of frequency bandwidth and likely forms of 
pilot compensation. 

Flight data instrumentation requirements are a direct function of what 
piloting tasks are to be considered. For each identifiable task or 
"outer" control loop the following data are necessary — either from direct 
measurement or by suitable estimation: 
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• The command loop state variable for the task 

• Its first and, If possible, second derivatives with 
respect to time 

• The primary "control" 

• Any states which may be associated with Intermediary or 
"Inner" control loops needed for performing the task. 
There may even be alternative competing candidates for 
inner loops . 


Simulator-generated data are more likely to be complete, accurate, and 
noise-free, but flight data will usually suffer omissions and distortions. 

The sampling rate requirements for NIPIP depend upon the bandwidth of 
the loop being examined. Solutions for outer- (task) loop pilot control 
strategy or task execution dynamics should normally require less frequent 
sampling than for inner-loop characteristics. Where inner- and outer-loop 
characteristics are estimated simultaneously (as in the previous approach 
cases) then the inner-loop bandwidth should dictate sampling rate. The 
rate of 50 saraples/sec was found adequate for successful analysis of heli- 
copter maneuvering 

Nonlinearities related to quantization or roundoff of recorded data 
should be viewed with concern. Double precision (e.g., coarse channel 
plus fine channel) may be necessary for any states crucial to a given 
piloting task. The quantization bands of 9 and <j> for the DFBW F-8 
(0.1 deg and 0.02 deg) might be used as guides for unacceptable and ac- 
ceptable coarseness of attitude angles, respectively. The 30 ft 
quantization of height was unacceptable for identification in the flare 
task and precluded analysis of the formation flight task. 

It is also recommended that careful consideration be given beforehand 
to data reconstruction and estimation schemes for any important state 
variables which cannot be directly measured and/or recorded, because such 
advance consideration may well exert an influence on the repertory of 
variables that can be measured and recorded . 


TR-1 188-2 


92 



B. AFTI/F-16 APPLICATIONS 


Simulation of the overall task, pilot, and vehicle Is an excellent way 
to verify the NIPIP outputs. That is, use the NIPIP outputs to simulate 
the pilot's control strategy and then compare the simulated outputs of the 
task, pilot, and vehicle to the actual outputs. 

One particularly attractive target for NIPIP application is the 
Advanced Fighter Technology Integration (AFT)l/F-16 flight program. Its 
concern with how to use the many varied mission-oriented flight control 
modes makes direct measurement of pilot control strategy and task execu- 
tion an appealing option. 

The following excerpt, taken from Ref. 17, provides a short background 
description of the AFTI/F-16 and Its program objectives. 


"The Advanced Fighter Technology (AFTI)/F-16 program is 
in response to today's European scenario, characterized by 
increased numbers of enemy targets both on the ground and 
in the air and an increasingly hostile air space surround- 
ing these targets . This changing environment required 
timely improvements in present IJSAF fighter lethality and 
survivability. The primary and continuing objective of 
the AFTI program, co-sponsored by the Air Force, NASA, and 
Navy is to provide for the development, integration, 
flight evaluation, and demonstration of emerging fighter 
technologies, and transition of the integrated technol- 
ogies to future system applications. The AFTI Fighter 
Attack Technology (AFTI/F-16) program will develop, inte- 
grate, and flight test a set of technologies to improve 
the survivability and weapon delivery accuracy of tactical 
fighters in air-to-air and air-to-ground attacks, through 
integration of advanced technologies into a single seat 
demonstrator vehicle which permits a realistic evaluation 
of technology benefits, penalties, and overall mission 
effectiveness . 

"The AFTI/F-16 vehicle has particular importance as a 
long life demonstrator aircraft with the flexibility, 
versatility, and capability in terms of performance and 
systems to serve as a future technology development test- 
bed. A full-scale development F-16 aircraft is the test 
vehicle. Extensive modifications were made for 
installation of a sophisticated data instrumentation 
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system, modified inlet with canards, new flight control 
system, and a dorsal fairing to accommodate the instrumen- 
tation equipment. Additional information on the AFTI/F-16 
can be found in Ref. [18]. 

"The overall objective of the AFTI/F-16 Advanced 
Development Program is to demonstrate separately, and in 
combination, advanced fighter technologies to improve air- 
to-air (AA) and air-to-surface (AS) weapon delivery 
accuracy and survivability. These technologies include a 
Digital Flight Control System (DFCS), Automated 
Maneuvering Attack System (AMAS) , pilot/vehicle interface 
(PVI) advancements, and advanced task-tailored control 
modes utilizing direct force control and weapon line 
pointing. Development, integration, and flight validation 
of these fighter attack technologies have been separated 
into DFCS and AMAS program phases. 

"The DFCS is a full-authority, triplex, digital fly- 
by-wire flight control system. The DFCS is mechanized to 
implement task-tailored manual control modes, including 
decoupled (six independent degrees of freedom of control- 
configured vehicle) flight control. Figure [49] shows 
that the pilot need only push a button to change the 
functions of cockpit controllers and displays. For the 
AMAS phase, the effective utilization of the advanced 
technologies requires the integration (coupling) of the 
fire and flight control functions. The integrated system 
will tie together a director fire control system, an 
advanced sensor-tracker, and the flight control system to 
provide precise automated weapon line control and weapons 
delivery. With the coupled system the azimuth and eleva- 
tion fuselage pointing capability of the aircraft provides 
an expanded envelope of fire control solutions; i.e., an 
enlarged pipper. The pilot need only capture the target 
within the expanded pipper envelope, and the fire control 
system will automatically command the flight control 
system to null aiming errors to assure a hit. This con- 
cept will profoundly influence fighter effectiveness in 
both AA and AS missions. 

"Pilot/vehicle interface advancements will be incor- 
porated to provide crew station capabilities and 
environment commensurate with the increase in total ve- 
hicle capabilities provided by the other technologies in 
each phase. The DFCS phase will focus on core technology 
development. The technologies of prime interest will be 
manual flight path control, avionics integration, and 
advanced controllers and displays. In the AMAS phase the 
allocation of function between the pilot and vehicle will 
be redistributed as a result of the DFCS experience. 
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Figure 4-9- AFTl/F-1 6 Multimode Flight Controller Commands 

(From Ref. 1 7) 
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Those tasks best performed by the machine will be auto- 
mated. Technological advances in sensors, fire control 
modes, and weapons fusing will be integrated with the DFCS 
capabilities . 

"An example of advanced technology integration and 
utilization is in the AMAS precision low altitude maneu- 
vering attack scenario . The technologies involved in this 
scenario include: 


1 . Flight path control with full authority 
digital flight control. 

2. Task automation with integrated flight and 
fire control and low altitude radar 
autopilot . 

3 . Advanced sensor-tracker with low drag FLIR 
and laser ranger installation. 

4. Integrated avionics and weapons fusing. 

5 . Cockpit development including multi-purpose 
displays, wide field of view heads up 
display, helmet-mounted sight and voice 
command . 

6. Weapons interface with pilot consent and 
auto-release . 


"These technologies together give the AFTl/F-16 the 
ability to more effectively attack ground targets. A low 
altitude radar autopilot allows survivable ingress and 
egress . AMAS automated air-to-surface bombing modes 
provide the capability for flexible target acquisition, 
precise tracking, automated ingress/attack steering, and 
automated weapon release for both low altitude, or stand- 
off delivery direct, or high-g turning attacks." 


Because of the AFTI program's emphasis on how a pilot uses the numer- 
ous flight control modes, task and pilot control strategy measurement 
offers a useful kind of documentation. There is the potential for detect- 
ing subtle differences in control strategy from one mode to another which 
could signal display deficiencies, natural pilot-to-pilot or run-to-run 
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variations, relative success in task execution dynamics, and relative 
distribution of pilot workload among task components. 

In order to succeed in pilot identification, however, the foregoing 
analysis cases point up the requirement for high quality flight data. 
This lesson should therefore play a key role in evaluating AFTI/F-16 
needs . 

Properly manipulated, NIPIP can be used for any of the basic tasks and 
maneuvers connected with the AFTI flight testing. The task and maneuver 
descriptions contained in Appendix A of Ref. 17 serve as a starting point 
for establishing command loops and primary flight controls. 

An example is shown below for the "air-to-surface tracking, bomb" 
maneuver defined in Appendix A of Ref. 17. 


Air-to-Surface Tracking, Bomb 

1. Set-up inbound to the target at 3500 ft above ground level. 

2. Upon reaching the point where the target is 10 deg below the 
horizon, pushover and track the target with the flight path 
marker. 

3. Use only the controllers specified in the run table. 

4. Recover from the dive at a safe altitude. 

The corresponding configuration and flight condition run table is pre- 
sented in Table 5. For Run SC-564, the decoupled bombing mode would be 
selected and stick and pedal controls used (direct lift control via throt- 
tle would not be available). Thus: 


Sp (stick, pitch axis) for flight path maneuver 
enhancement 

S R (stick, roll axis) for roll rate 
P (rudder pedals) for flat turn 
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TABLE 5 


CONFIGURATIONS AND FLIGHT CONDITIONS 


RUN NO. 

AS/ 

MACH 

ALT. 

CONFIG. 

EXT. 

LOAD 

DFCS 

MODE 

PRIORITY 


MANEUVER 



SC- 560 

400 

kcas 

3. 

5K 

C 

R 

A/ 

’S 

SASB 

] 

l 

Air-to-surface 

bomb 

tracking. 

stick 

only 

SC-561 









SASB 



Air-to-surface 

bomb 

tracking, 

stick 

and pedals 

SC-562 









DASB 



Air-to-surface 

bomb 

tracking, 

stick 

only 

SC-563 


r 







DASB 



Air-to-surface 

bomb 

tracking. 

stick 

and throttle 

SC-564 

400 

kcas 







DASB 



Air-to-surface 

bomb 

tracking, 

stick 

and pedals 

SC-565 

500 

kcas 







SASB 



Air-to-surface 

bomb 

tracking, 

stick 

only 

SC-566 









SASB 



Air-to-surface 

bomb 

tracking, 

stick 

and pedals 

SC-567 









DASB 



Air-to-surface 

bomb 

tracking, 

stick 

only 

SC-568 

\ 

r 

1 

f 


f 

1 

r 

DASB 

1 

r 

Air-to-surface 

bomb 

tracking, 

stick 

and throttle 

SC-569 

500 

kcas 

3. 

5K 

C 

:r 

A/S 

DASB 

• 

1 

Air-to-surface 

bomb 

tracking. 

sitck 

and pedals 



Task segments implied are: 


1. Inbound to target, level at 3500 ft above ground level 

2. Pushover and track target with flight path marker (HUD 
symbol) 

3. Recover from dive. 


For each segment, the implied command loop/control combination is: 


1 . 


2 . 


3. 


Inbound , 

level 

flight 


h > 

0 c i 

6 + S p 

(h 

y - 

V 

*c 5 

* ♦ S R 

(y 

Pushover, 

track 

target 



Recovery 

9 + S p 

*c + S R 


A 

= height) 

A 

= lateral path displacement) 


NIPIP would therefore require definition of a finite difference equation 
for each task or control strategy structure implied by the above 
combinations . 

For example, for h -*■ Sp the closed-loop task dynamics might reasonably 
be given by a second-order characteristic equation: 


h + 2C^h + w h = 0 


(35) 
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Hence the finite difference equation would be: 


h(n) = -2?toh(n) - w h(n) + bias 


( 36 ) 


where h(n) and h(n) can be estimated from a z and h as described in 
Appendix A. 


Solving for 2C“, w , and the bias provides an estimate of closed-loop 
activity in holding altitude. 


w aggressiveness in altitude regulation 
C damping, freedom from PIO 


Carrying this example further, pilot control strategy in the same 
altitude loop could be measured by considering the correlation between the 
control Sp and the inner- and outer-loop states 6 and h. The same dif- 
ferential equation form demonstrated in the previous examples (Cases 2 and 
5) would be appropriate. 

Note that only one of the bombing segments has an inner- and outer- 
loop combination. The tracking and recovery segments probably involve 
only inner loops. Nevertheless there are features worthy of study. For 
example, what does the pilot do with the lateral stick during the tracking 
segment? Is there stick and pedal coordination? Is such coordination 
subliminal or does the pilot consciously apply it? 
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Suggested candidates for NTPIP difference equations for the target 
tracking and recovery segments are: 


Target Tracking, Vertical Axis: 

S p (n) = k 1 Sp(n-l) + k 2 e y (n) + k 3 e v ( n -l) + k^ 
Target Tracking, Lateral Axis: 

P(n) = k^P(n-l) + k 2 ^(n) + k^e^Cn-l) + k^ 
and, assuming some coordination with lateral stick, 
S R (n) = k 3 S R (n-l) + k^P(n) + k^ 

Recovery, Vertical Axis: 

Sp(n) = k^SpCn-l) + k 2 9(n) + k 3 0(n~l) + k^ 
Recovery, Lateral Axis: 

S R (n) = kjSpCn-l) + k 2 <Kn) + k 3 <|>(n-l) + k^ 


These forms provide for identification of pilot lead and lag (or delay) 
compensation along with general loop tightness. The difference equation 
forms can be altered to enhance the definition of any of these specific 
qualities where desired. For example, additional degrees of freedom 
involving the second (or more) previous sample(s) in the "controller" 
terms will better define lag characteristics. 

In other instances, if anticipated loop bandwidth permits, the analyst 
may incorporate only the second (or mth, where m is an integer) previous 
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saraple(s) in the "controller" terms to improve definition of lag 
characteristics . 

Pilot questionnaires and briefing procedures should be designed to aid 
in the task and control strategy identification process. At the same time 
limitations in the pilot's ability to analyze control strategy or task 
execution introspectively should be appreciated. 

The main factors to probe in connection with any task are the choice 
of controls, how tasks are segmented, and what cues are used. These ques- 
tions may be aided by helping the pilot subject to construct conventional 
control loop block diagrams. It may also be instructive to the analyst to 
ask the pilot about special "tricks" in his control strategy such as co- 
ordination of two controls, anticipation, or use of unusual kinds of 
cues. Finally, it is important to determine any factors which might tend 
to make a given run atypical. 


C. AUTOMATIC SELECTION OF PILOT CONTROL STRATEGIES 

Provisions for automatic pilot control strategy identification were 
implemented in the version of NIPIP documented in Ref. 7. These consisted 
of multiple simultaneous pilot control strategy difference equation solu- 
tions along with conventional goodness of fit metrics. This permits on- 
line assessment of NIPIP results in either a flight or simulation 
environment. Comparisons can be made in terms of several parameters de- 
pending upon how the analyst chooses to specify the NIPIP difference 
equation options. 

It must be stressed, however, that truly "automatic" pilot control 
strategy selection is fraught with hazards and unknown consequences at 
this stage. Control strategy selection must really be accomplished in a 
manual, interactive mode using engineering judgment and the results of 
past experience. With this strong caveat, we shall now expand on how the 
limited "automatic" selection tools might be exploited. 
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There are essentially two stages of pilot control strategy identifica- 
tion where the above NIPIP features can be effectively used. One is 
connected with basic control loop structure identification, the other with 
selection of control compensation identification or data smoothing forms, 
examples of which were given for the target tracking task in the previous 
topic. The order of these two steps is not clear — both may be done at 
once, in fact. 

For control loop structure identification, several competing NIPIP 
difference equations might be chosen using different combinations of pri- 
mary controls and feedback variables. For example: 


"Frontside" : 


5 e = 


k 0 9 + k£h + 


k h h 


+ *>; 


6 T “ 


k u u 


+ c 


"Backside" : 


6 e = k e 0 + k^u + kffu + b; 6 T = + k^h + c 


"Backside" with ^ 

throttle coordination: e 


k 0 9 + k u u + k x 6 T + b; 6 T = k^h + k^h + c 


It is thus possible to distinguish the best choice of control structure by 
observing any of the available goodness-of-f it indicators either mentioned 
previously or any of those which will be suggested in the computer- 
graphics discussion. 


D. INTERACTIVE COMPUTER GRAPHICS 

Because, at this stage, pilot control strategy identification is an 
iterative process, it is desirable to have the means for quick, effective 
evaluation of NIPIP results. The version of NIPIP now operational 
produces a large array of tabulated calculations, but these require a 
separate processing in order to fully interpret their quality. 
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The use of an interactive computer graphics scheme directly tied to 
NIPIP would be helpful indeed. Some experience in this area has been 
gained in previous programs, and much was learned in the study reported 
here to serve as a basis for recommendations. The following paragraphs 
describe these recommendations. 

In general, the dynamic response of systems can be presented in many 
forms, each providing its own special insights. This can include the 
domains of time, frequency and phase plane, continuous or discrete. None 
of these alone can be regarded as wholly adequate for the analyst, 

though. It is advisable to exploit as many separate presentations as 
possible for the purposes of finding an acceptable solution and for con- 
firming it, 

This report presents some of the ways of portraying NIPIP results, but 

it is a fairly limited sample. The recommendations of this section con- 

tain many more possibilities even though not all have been tested for 
their effectiveness. 

Interactive computer graphics, to be effective in the NIPIP role, must 
be sufficiently flexible to accommodate several kinds of presentations, 

reasonably high resolution, fast enough to keep up with a running NIPIP 
solution (which could be on-line, real-time) and able to generate a hard 
copy if desired by the analyst. Each of these attributes will be dis- 
cussed, in turn, in the following paragraphs. 

First, the NIPIP user Is concerned with observing (a) the data being 
analyzed and (b) the solution in its various alternative forms. The for- 
mer provides a starting point for assuming a candidate loop structure 
form, the latter, the adequacy of the solution and insight for refine- 

ments. Hence a computer graphics scheme needs direct access to both the 
input to and the output from NIPIP. 

The flexibility required in plotting relates to choice of independent 

and dependent variables and to scaling. There can be no hard and fast 

rules. For Inner-loop concerns, time scales might be expanded and choice 
of state variable limited to inner-loop quantities — pitch attitude, roll 
attitude, yaw, heading, and sometimes vertical velocity or flight path 
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angle. Outer loops would necessitate another set of plots. While some 
specific plotting objectives will be given shortly, there should always be 
the ability to modify them. 

The form of computer graphics most useful to the NIPIP user is a hard 
copy, scaled, two-dimensional plot. It would be convenient, however, to 
use a CRT display as an intermediate step in obtaining a hard copy. 

Several basic plots of input data and NIPIP solutions are presented in 
Fig. 50. Each is discussed below: 


1. Control and State Variable Time History 

As a first step in the pilot identification process, it is useful to 
inspect simple time histories of the command loop state variable and the 
suspected control for that state. Other states and controls may also be 
of interest, however. Further it is beneficial to superimpose these time 
histories in order to gain insight about correlation, phasing, relative 
frequency content, and task segmentation. 

A computer-graphics display of raw data time histories may require 
positive labeling of individual states. This could be difficult to ac- 
complish via conventional line coding (solid, broken, or alphanumeric 
symbols). A multicolor display would be feasible, however, for both a 
hard copy plotter and a CRT. A variety of multicolor plotters are on the 
market at reasonable cost (e.g., Huston Instruments and Hewlett- 
Packard). Color monitors are also available and easily driven by low cost 
microcomputers such as Apple or TRS-80. The main difficulty in using 
color media lies in the cost of reproduction of large numbers for desem- 
ination of reports. While color xerography is readily available, it is 
expensive . 
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Assumed Control Strategy: 



Basic Finite Difference Equation Used in NIPIP : 




Figure 50 . Recommendations for Interactive Graphics with NIPIP 
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k. Phase Plane Comparison 



5 . Time -Frequency Describing Functions 




Figure 50 (Concluded) 
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2. Control and State Variable Phase Plane 

The same data plotted as time histories can also be plotted without 
time as the independent variable. The value of phase planes is that cor- 
relation between pairs of variables is easily seen, and non-linearities 
can be detected and even identified. (Ref. 15 contains a large catalog of 
phase planes for nonlinear elements.) 


3. Time History Comparison of 

Model Reconstruction with Raw Data 

One rather clear way of judging goodness-of-model is to reconstruct a 
control or state using a set of model coefficients obtained by NIPIP. For 

A 

example, if NIPIP solved for the coefficients 3 and b assuming the dif- 
ference equation: 


6 = a 0 + b 

n n 


(37) 


then a modeled "6" defined as 


6 = a 9 + b 

n n 


(38) 


can be generated using the raw data for 0. This <5, in turn, can be com- 
pared directly with the actual 5 in order to help to confirm the model. 
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4. Phase Plane Comparison of 

Model Reconstruction with Raw Data 


This is a counterpart to the previous graphic form where <5 and 6 are 
plotted . 


5. Time-Varying Describing Functions 

One interesting visualization of NIPIP results is the construction of 
a time history of the frequency response, i.e., gain and phase as in 
Section IV of this report. This concept lacks mathematical rigor, but it 
does help to evaluate the consistency and general character of a NIPIP 
solution. 
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APPENDIX A 


SINK RATE ESTIMATION 


Since the flight test data did not contain any direct measurement of 
approach slope or sink rate, it is necessary to provide a reasonable 
estimate of this parameter. Thus, as a first step in the analysis, it is 
necessary to use the existing flight data to estimate sink rate and to 
determine the sink rate command profile for the maneuver. A constraint in 
the choice of methods was that this report was to define the NIPIP 
analysis procedure, not techniques in state variable estimation. 

Complimentary filtering was used in estimating sink rate to take 
advantage of the data available. The altitude data is appropriate for 
low-frequency estimation of sink rate while vertical acceleration is 
appropriate for high frequencies. Complimentary filtering allows the data 
to be combined in a way that takes advantage of these relative 
strengths. The complimentary filter in continuous form is: 


h 


— §5-h +-J- 

s + a s + a 


h 


(37) 


where 

s is the Laplace operator 

h is the measured altitude 

h is vertical acceleration estimated from 
measured normal acceleration 

h is the estimated sink rate 

and a is the characteristic frequency of the filter 
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The complimentary filter was implemented in finite difference equation 
form as: 


h 

n 



+ ah 

n 


ah 


n-1 




(38) 


The characteristic frequency, a, was determined empirically to accommodate 
the sample period as well as the quantization in the measured altitude and 
the noise in the vertical acceleration. A value of 0.1 was chosen for 
use. Larger values produced a sink rate estimate which showed the 
quantization of the measured altitude. 

The vertical acceleration used in the complimentary filter was 
estimated from the flight data using two separate approaches. The first 
approach used all of the measured aircraft states, both lateral and 
longitudinal, to reconstruct the vertical acceleration. Assuming the 
normal acceleration to be measured at the center of gravity, one can 
estimate the vertical acceleration by 


h = -(-a^ + g cos 9 cos <j> - PV + QU) cos 0 cos <f> (39) 


where a n is the measured normal acceleration. The second approach used a 
simplified method which corrected measured normal acceleration only for 
pitch attitude effects to obtain vertical acceleration. Vertical 
acceleration in this case is given by 

h = [(a - g cos 9 1 + g sin 00+ QU 1 cos 9 (40) 

lv n ° o' o <J o 

where cos 0 Q was taken as unity. 
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This latter estimation approach is appropriate when large lateral 
maneuvers are not present. However, as there are large lateral maneuvers 
in the flight data (as seen in the flight data time histories) the first 
method was also relied upon to estimate the vertical acceleration. 

The final aircraft state to be estimated is h . This state was 
estimated by using pitch rate and vertical acceleration by: 


h 



(41) 


where T fl was determined to be 
°2 


equal to 1 sec. 


At this point something should also be said about the quality of the 
data which is to be used in the estimation technique. It goes without 
saying that the better the data the better the chance of success and the 
more reliable the outcome should be. However, there is a point at which 
the identification technique cannot be relied upon to provide an accurate 
estimate due to poor quality input data. This aspect of the quality of 
the input data will be addressed in the numerical results in both 
Section IV and Appendix B of the text. 
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APPENDIX B 


INVESTIGATION OF THE EFFECTS OF QUANTIZATION IN PITCH ATTITUDE 
ON THE IDENTIFICATION OF PITCH ATTITUDE AND SINK RATE CONTROL STRATEGY 


A simplified simulation of a generic F-8 aircraft was used to explore 
pitch attitude quantization problems encountered with the DFRC flight 
data. The features of the most concern were the approximately 0.1 deg 
incremental steps in the pitch attitude data channel. This STI F-8 simu- 
lation involved defining an autopilot which would allow for changes in the 
flight task during the computation of a time history. The autopilot de- 
rived for this study commanded a steady rate of climb or decent. A block 
diagram of the generic F-8 aircraft simulation and autopilot is given in 
Fig. 51. This particular autopilot was designed to mimic the pilot's 
technique for controlling an approach to landing, which is one of the 
flight tasks provided by DFRF in Section IV of this report. 


The inner-loop pilot strategy, Y , and the outer-loop pilot strategy, 

^ 0 

Yp^, were chosen to be of the form Kq and K^/s, respectively, where K^/s 
refers to the Laplace transform of an integrator. The gains, Kq and K^, 
are equal to 1.675 rad/rad and 0.0005 rad/ft, respectively, which gave 
inner- and outer-loop crossover frequencies of 3.0 and 0.2 rad/sec, 
respectively. 


The change in flight task was demonstrated by commanding the aircraft 
from an altitude hold (i.e., straight and level flight) to a steady rate 
of descent (i.e., pushover or flight-path hold) flight task as shown in 
Fig. 52. The quantization of pitch attitude in 0.1 deg incremental steps 
is also shown in Fig. 52. 

Figure 53 presents the results of using a longer sampling period and 
starting the identification procedure at the initiation of the push-over 
(T = 15 sec). Three cases are shown: (a) a sample period of 0.1 sec with 
0.1 deg quantization, (b) a sample period of 0.1 sec with no quantization, 
and (c) a sample period of 1 sec with 0.1 deg quantization. All three 
cases have a time window of 25 sec and use five degrees of freedom for 
estimation (see Eq . 10 in Section IV). As seen in the 0.1 sec sampling 
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Bandwidth of 6 -loop is 3.0 

Bandwidth of h-loop is 0.1 

. o.o for 0 < t 1 p sec 

h - 

-d'7-5 fps for t _> 1 3 s 

Figure 51 • 











Figure 52 . Pushover Time History of Pitch Attitude 
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period traces where quantization is present, merely starting the window 

after the push-over does not improve the ability to identify the control- 

loop elements. However starting after the push-over and increasing the 

sample period to 1 sec improves the identified solution. The outer-loop 

identified solution, Y . , with a 1 sec sample period and a 0.1 deg quanti- 

"h 

zation in pitch attitude essentially matches the actual solution 

identified with a 0.1 sec sample period and no quantization. The inner- 

loop identified solution, Y , for a sample period of 1 sec with 0.1 deg 

P0 

quantization does not match the actual inner-loop solution with a sample 
period of 0.1 sec and no quantization, because the Nyquist frequency is 
approximately equal to the crossover frequency of the inner loop. 


Increasing the sampling period from 0.1 to 1 sec does improve the 

identified solution for Y . as seen by comparing the 0.1 sec and 1 sec 

' J h 

sampling period traces with 0.1 quantization in pitch attitude. It should 
be noted that usually the Nyquist frequency must be greater than the 
crossover frequency to obtain accurate estimates; however, in this case 

the effects of quantization degrade the ability of the procedure to iden- 
tify the inner-loop control strategy. Thus using still longer sampling 

periods does not improve the ability to identify the inner-loop but does 
at the same time improve the ability to identify the outer-loop control 
strategy. Shorter sampling periods as shown by the 0.1 sec sample period 
degrade the ability to identify both the inner and outer control-loop 
elements in the presence of the specified quantization. Hence it is not 
possible to identify accurately the inner control-loop strategy in the 
presence of this particular level of quantization. 


The preceeding results show the effects of quantization in an inner 
control-loop variable on both inner- and outer-loop strategy identifica- 
tion. The effects of pitch attitude quantization were shown to degrade 
the ability to identify the control-loop elements. However it was also 
shown that when the quantization is present only in the inner-loop, it was 
still possible to identify the outer-loop control strategy by adopting a 
longer sampling interval. The results of this investigation support the 
initial conclusion that quantization of pitch attitude does, in fact, 
degrade the ability of the NIPIP to correctly identify the inner-loop 
control strategy. 
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